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Using variational perturbation theory, we calculate the shift in the critical temperature Tc up to 
five loops to lowest order in the scattering length a and find ATc/Tc'^'' ~ (1.14 ± 0.11) an^^^, where 
'^f^ • n is the particle density. Our result is lower than the latest Monte Carlo result (1.32 ± 0.02) an^^^. 
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bJ)' The effect of a small repulsive interaction upon the critical temperature Tc of a Bose-Einstein condensate 
^ ' (BEC) has been a matter of controversy for many years, and the various results are converging only very 

slowly towards a common answer. Here we want to contribute to the ongoing discussion with a further result 

obtained from a resummation of Feynman diagrams up to five loops. 
. The interacting Bose-Einstein condensate is described by the euclidean action 
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^^ = J^ dry d^x|r(x,r) (^a,-^V^jV'(x,r)-/iV'*(x,r)^(x,T) + :^[^(x,T)^*(x,r)]'| , (1) 

, where M is the mass of the bosons, /3 the inverse temperature in natural units with h = ks — 1, a is the 
s-wave scattering length, and /Lt the chemical potential. The free system has a transition temperature 
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^ ' where n is the particle density. A small relative shift of Tc with respect to ri"'' can be calculated from the 
general formula 
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where n^^^ is the particle density in the free condensate and An its change at Tc caused by the small 
interaction. For small a, this behaves like P, Q 
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where 4 = -647rC(l/2)/3C(3/2)5/3 ~ 19.7518 can be calculated perturbatively, whereas c\ and C2 require 
nonperturbative techniques since infrared divergences at Tc make them basically strong-coupling results. 
The standard technique to reach this regime is based on a resummation of perturbation expansions using 
the renormalization group first applied in this context by Ref. Jsj. Recently, however, it has been 

shown by calculating the best known critical exponent a of superfluid helium from Satellite experiments 
that the accuracy of strong-coupling results can be surpassed by much simpler variational perturbation 
theory Mi. 

Up to now, C2 has been inferred only from Monte Carlo data to be C2 ~ 75.7 ± 0.4. In order to find the 
leading coefficient c\ , one may take advantage of an important simplification due to the fact that An can be 
calculated from the classical limit of the field theory, which is governed by the three-dimensional action 

As^^pJ d'x |^o*(x) (-^V^ - f?j V'o(x) + ^ [^„*(x)V'o(x)]'| . (5) 

This is a special case = 2 of the more general 0(A)-invariant 0** field theory 
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where the A^'-component field ({> = (f>2, ■ ■ ■ , 4>n) is related to the original field -0 for A'' = 2 by tpix) = 
VmT[0i(x) +i<^2(x)]. The square mass is = — 2M/z, and the quartic coupling is w = ASnaMT. Using 
this relation, the shift of the critical temperature © can be found from the formula 
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corresponding in Eq. (0J) to 
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The three-dimensional theory is superrenormalizable and requires only mass counterterms which shift the 
original bare mass m to the renormalized mass m^. A calculation of the Feynman diagrams in Fig. ^yields 
the following five- loop perturbation expansion for the expectation value {(fP' /u) |10lllll| 
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where 02 = log(4/3)/2 w 0.143841 and the other constants are only known numerically [l^ : 

03 = 0.642144, 041 = -0.115069, 042 = 3.128107, a43 = 1.63, 044 = -0.624638, 045 = 2.39. 
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Writing the above expansion up to the Lth term as Fl{u) ~ I]j^_j^//(u/47rTO,.)', the expansion coefficients 
for the relevant number of components N = 2 are 12]: 
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-126.651 X 10"^ /o = 0, A = -4.04837 X 10"^ /a = 2.39701 x 10"^ /g = -1.80 x 10"^ (11) 



We need the value of the series Fl(u) in the critical limit rrir — *■ 0, which is obviously equivalent to the 
strong-coupling limit of Fl{u). As mentioned above, this limit should be most accurately found with the 
help of variational perturbation theory 0, 0, 13 • 

If the series were of quantum mechanical origin, we could find this limit by applying the rules of naive 
variational perturbation theory We form the sequence of truncated expansions Fl{u) for 1,2,3 and 

replace each term 

(^./m.)'->X'[l-l]-i/f (12) 
where the symbol [1 — 1]^ is defined as the binomial expansion of (1 — l)*" truncated after the fcth term 
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FIG. 1: Diagrams contributing to the expectation value ((^^). 



3 



TABLE I; Trial functions for the naive quantum-mechanical variational perturbation expansion 
y^QM ^ _o.05968317i'-i - 0.0000322159 if, 
W^^ = -0.04973597^"^ - 0.0000483239 K + 1.51792 10"^ , 
W^^ = -0.04351897^-^ - 0.0000604049 K + 3.03584 IQ-^ - .908 10"^ K^. 



The resulting expressions must be optimized in the variational parameter K. They are listed in Table U 
The approximants W^^^ have extrema W^^^""^ « -0.00277, +0.00405, -0.0029, corresponding, via ©, to 
ci « 3.059, —4.46, 3.01. These values have previously been obtained in Ref. in a much more complicated 
way via a so-called ^-expansion. Note the negative sign of the second approximation arising from the fact 
that an extremum exists only at negative K. According to our rules of variational perturbation theory one 
should, in this case, use the saddle point at positive K which would yield = —0.00153 corresponding to 

ci ~ 1.69 rather than -4.46, leading to the more reasonable approximation sequence ci « 3.059, 1.69, 3.01, 
which shows no sign of convergence. In 1^3^"^, there is also a pair of complex extrema from which the authors 
of Ref. ^3 extract the real part Re W3^j,^^jjipi(,x ~ —0.00134 corresponding to ci « 1.48, which they state as 
their final result. There is, however, no acceptable theoretical justification for such a choice [l^ . 

This lack of convergence is not astonishing since naive quantum mechanical variational perturbation theory 
is inapplicable to field theory, contrary to ubiquitous statements in the literature 15]. A simple but essential 
modification is necessary to allow for the well-known fact that there are anomalous dimensions in the critical 
regime of fluctuating fields. This modification was discovered in Ref. and tested by the fact it reproduces 
in D = 4 — e dimensions exactly the known e expansions of renormalization group theory j^. In Z? = 3 
dimensions, it leads to the most accurate critical exponents so far (see in particular Chapters 19 and 20 in 
the textbook 

The correct procedure goes as follows: We form the logarithmic derivative of the expansion (Q: 



OlogU /_i \mrj f-1 \mrj V /-I /-l/ 



In order for F{u) to go to a constant in the critical limit rrir — » 0, this function must go to zero in the 
strong-coupling limit u — s- 00. Writing the expansion as Pl (w) = —1 + S/i^2 hiu / 4:TTmry , the coefficients are 

&2 = 0.0639293, 63 = -0.056778, 64 = 0.0548799. (15) 

The sums Pl(u) have to be evaluated for u — s- cx) allowing for the universal anomalous dimension u by which 
the physical observables of (/>^-theories approach the scaling limit 0, Q ■ The approach to the critical point 
A + B{mr/u)'^ where uj' — — 77/2) 16]. The exponent 77 is the small anomalous dimension of the field 
while Lo is related to the famous Wegner exponent 5] of renormalization group theory A = ujv. Here it 
appears in the variational expression for the strong-coupling limit which is found 0,0] by replacing (u/m^)' 

by K^l - where q = 2/uj'. Thus we obtain the variational expressions 

= -1+(IIi + ^)k' + ^K' (16) 
2/1 , 3/ig , /3/2 , 9/29\ ,^3 , f-2f,^ , 4/3 



The first has a vanishing extremum at u'^ = 0.592, the second has neither an extremum nor a saddle point. 
However, a complex pair of extrema lies reasonably close to the real axis at us'^ = 0.635 ± 0.116, whose real 
part is not far from the true exponent of approach uj'^ « 0.81 0,0], to which cj^ will converge for order 
L ^ 00 0. Given these w'-values, we now form the variational expressions Wl from by the replacement 
{u/mrY — > if']l — i]]^'!'^!^, which are 

W2 = f-i (1 - ^9 + J?') + hK, (18) 
- /-I (1 - ^9 + - ^9') + (1 + I) ^ + ^2^'' (19) 
^ /-i(l-i'Z + |.^-4.^+3^.^)if- + /t(l + J. + i.^)i^ + /.(l + .)i^^ + /3i^^ (20) 
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The lowest function W2 is optimized with the naive growth parameter q = 1 since to this order no 
anomalous value can be determined from the zero of the beta function (|14|) . The optimal result is W2^^ = 
— •\/log[4/3]/6/87r^ « —0.00277 corresponding to ci = 3.06. The next function W3 is optimized with the 
above determined 93 = and yields W!^''^^ ~ —0.000976 corresponding to ci = 1.078. Although uj'^ is not 
real we shall insert its real part into W4 and find W^^^ = —0.000957 corresponding to ci = 1.057. The three 
values of ci for L = L — 1 = 1,2,3 can well be fitted by a function ci w 1.053 + (see Fig. Such a 

fit is suggested by the general large-L behavior a + her'^^^ ^ which was derived in Refs. Due to the 

smallness of 1 — ~ 0.2, this can be replaced by sa a' + b' /L^ . 

Alternatively, we may optimize the functions VFi^2.3 using the known precise value of q^c — w 
2/0.81. Then W2 turns out to have no optimum, whereas the others yield W^'^^ « -0.000554, -0.000735, 
corresponding via Eq. Q to ci = 0.580, 0.773. If these two values are fitted by the same inverse power of 
L, we find ci ~ 0.83 — 14/L^. From the extrapolations to infinite order we estimate ci^oo ~ 0-92 ± 0.13. 




0-5 ..•••'■*ci Si 0.830 - 

1 1.5 2 2.5 3 3.5 ^ 

FIG. 2: The three approximants for ci plotted against the order of variational approximation L = L — 1 = 1, 2, 3, 
and extrapolation to the infinite-order limit. 



This result is to be compared with latest Monte Carlo data which estimate ci k, 1.32 ± 0.02 0, H^ . 
Previous theoretical estimates are ci « 2.90 2.33 from a l/jV-expansio n [20ll ). 1.71 from a next-to- 
leading order in a 1/A^-expansion j2j|, 3.059 from an inapplicable (5-expansion |22| to three loops, and 1.48 
from the same (5-expansion to five loops, with a questionable evaluation at a complex extremum and 
some wrong expansion coefficients (see [131 )■ Remarkably, our result lies close to the average between the 
latest and the first Monte Carlo result ci « 0.34 ± 0.03 in Ref. p^ . 

As a cross check of the reliability of our theory consider the result in the limit N 00. Here we must 
drop the first term in the expansion ^ which vanishes at the critical point (but would diverge for N 00 
at finite m^). The remaining expansion coefficients of (^(jP /uj /N in powers of Nu/Annir are 

A = -6.35917 10"^ A = 4.7315 10""^, A = -3.84146 lO"'^. (21) 

Using the N ^ 00 limit of uj' which is equal to 1 implying q = 2 in Eqs. (|19|l and H20(l . we obtain the two 
variational approximations 

= -0.00127183i\: + 0.00047315iv:^ = -0.00190775ii: + 0.00141945^^ - 0.000384146ii'^ (22) 

whose optima yield the approximations ci ~ 1.886 and 2.017, converging rapidly towards the exact large- iV 
resuh 2.33 of Ref. :20|, with a 10% error. 

Numerically, the first two 1 /A^-corrections found from a fit to large- Af results obtained by using the known 
large- A^ expression for w' = 1 - 8{8/3tt^N) + 2(104/3 - 9tt^ /2){8/3n'^N)'^ ^ produce a finite-TV correction 
factor (1 - 3.1/N + 30.3/A^2 + . . to be compared with (1 - 0.527/7V + . . .) obtained in Ref. 

Since the large- A^ results can only be obtained so well without the use of the first term we repeat the 
evaluations of the series at the physical value N — 2 without the first term, where the variational expressions 
for / are 

W2 - A (i + |)i^ + A^', 

W3 - A (1 + ^9 + K + f2{l + q)K^ + hK\ (23) 

The lowest order optimum lies now at = — /f (2 + (7)^/16/|, yielding ci = 0.942 for the exact q — 2/0.81. 

To next order, an optimal turning point of yields ci « 1.038. 
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At this order, we can derive a variational expression for the determination of lo' using the analog of Eq. (|14|l 
which reads 



_ d\ogF{u) ^^^h^^f^h_.fl 
dlogu fi rrir \ fi ff 

After the replacement (|12|l we find 

^ + + (25) 

/I V ji Ji 

whose vanishing extremum determines to' = 2/q as being 



2^2hh/ fi - 1 - 1 j ~ 0.675, (26) 

leading to ci w 1.238 from an optimal turning point of W^. There are now too few points to perform 
an extrapolation to infinite order. From the average of the two highest-order results we obtain our final 
estimate: ci « 1.14 ± 0.11, such that the critical temperature shift is 

AT 

^« (1.14 ± 0.11) (27) 
This lies reasonably close to the Monte Carlo number ci w 1.32 ± 0.02. 
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Note added in proof: 

Since this paper appeared on the Los Alamos server last October, the Feynman diagrams where recalculated 
more accurately and extended to six loops in Ref. 25] . The more accurate /a has increased by about 6% to 




0.5 ^,'*'' ci Si 0.914 - 2.7/L3 



FIG. 3: The four approximants for ci plotted against the order of variational approximation L = L — 1 = 1, 2, 3, 4, 
and extrapolation to the infinite-order limit. The second- highest aproximations are ci = 1.115 and 0.797 for self- 
consistent uj' = 0.745 and exact uj' ~ 0.81, the highest aproximations are ci — 0.959 and 0.886 for self-consistent 
uj' = 0.745 and exact uj' ~ 0.81, respectively. The latter dots cannot be distinguished. 

— 1.92 X lO^"', and the new expansion coefficient is /4 = 1.873 x 10^''. Using the new numbers to perform a 
variational evaluation fits of the first type including the anomalous first term /_i we find the curves shown 
in Fig. 13 yielding ci « 0.98 ± 0.06, only slighly larger than than the previous five-loop number 0.92 ± 0.13. 
The upper curve is now associated with the self-consistent w'-values — 0.626 and W4 = 0.745, the last 
being almost equal to the exact value 0.81. 

The evaluation of the second type where the anomalous term is omitted yields now ci = 0.942, 1.056, 1.120, 
for uj' = 0.81, and ci — 1.409, 1.384 for the variationally determined W3 = 0.621 and 0.638, as shown in in 
Ref. Combining the last two results one obtains the larger value ci = 1.23 ± 0.12, lying closer to the 

Monte Carlo number ci « 1.32 ± 0.02 than our five-loop result l|77|l . 
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